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We study dynamics of a mobile impurity moving in a one-dimensional quantum liquid. Such 
an impurity induces a strong non-linear depletion of the liquid around it. The dispersion relation 
of the combined object, called depleton, is a periodic function of its momentum with the period 
27m, where n is the mean density of the liquid. In the adiabatic approximation a constant external 
force acting on the impurity leads to the Bloch oscillations of the impurity around a fixed position. 
Dynamically such oscillations are accompanied by the radiation of energy in the form of phonons. 
The ensuing energy loss results in the uniform drift of the oscillation center. We derive exact results 
for the radiation-induced mobility as well as the thermal friction force in terms of the equilibrium 
dispersion relation of the dressed impurity (depleton). These results show that there is a wide range 
of external forces where the (drifted) Bloch oscillations exist and may be observed experimentally. 
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I. INTRODUCTION 

Recent experiments pji4j achieved a substantial progress in fabricating and studying dilute impurities immersed in 
one-dimensional (Id) quantum liquids. Such liquids are formed by ultracold bosonic or fermionic atomic gases placed 
in Id optical lattices [SHS] or using the magnetic confinement on atom chips [4] . A hyperfine state of a few atoms is 
then switched locally with the help of the RF magnetic field pulse [Hll]. As a result, the atoms become distinguishable 
from the rest of the quantum liquid and thus may be considered as mobile impurities. Another promising realization 
is achieved by placing ions of Yb^ , Ba^ or Rb"*" in the Bose- Einstein condensate of neutral *^Rb atoms [2 [3] . 

Remarkably, the impurities may be selectively acted upon with the help of external forces. In the case of neutral 
atoms the external force is gravity, uncompensated by the force from the vertical magnetic trap as the impurity 
atoms are created in magnetically untrapped states. In the case of an ion, the selective force is due to the applied 
electric field. These setups thus allow to drive the non-equilibrium dynamics of mobile impurities immersed in a Id 
quantum liquid. By releasing the trap after a certain delay time allows monitoring the resulting density and velocity 
distribution of impurities. 

It is important to emphasize that in this publication we are dealing with the case of impurities with finite bare mass 
M, not very different from the mass m of the particles in the background. The velocity of the impurity is therefore 
free to change as a result of absorbed momentum. This situation is qualitatively different from that of static, infinitely 
massive impurities [TMTU] . 

Dynamics of mobile impurities in quantum liquids is an old subject pioneered by Landau and Khalatnikov '10' and 
Bardeen, Baym, Pines and Ebner [Ulll^ in their studies of '^He atoms in superfluid ''He. These authors realized that 
at a finite temperature the impurities experience a viscous force from the normal component of the liquid, even if 
their speed is less than the critical superfluid velocity. The mechanism leading to such a viscous force was shown to 
be the Raman scattering of the thermal excitations of the superfluid i. e. thermal phonons scattering off the impurity 
■^He atoms. It was found [TOl [12] that the corresponding friction coefficient scales with temperature as in the low 
temperature limit. It was later discussed by Castro-Neto and Fisher [13j that in Id the Raman mechanism leads 
to the friction force proportional to due to the phase-space considerations. The same mechanism governs the 
temperature- induced acceleration of grey solitons [44] . Recently, two of the present authors [Ml [43] showed that in 
addition to phase-space effects the friction is extremely sensitive to the details of the interactions between impurity 
and the background liquid and vanishes for exactly solvable Id models. The friction therefore may be rather small in 
the spin-flip setup if the parameters are near the Yang-Gaudin [T3| integrable case. 

Due to the one-dimensional kinematics the slowly moving impurity tends to deplete the host liquid creating density 
and current gradients in its vicinity. At small momenta the dispersion of the impurity can still be described by 
quadratic law E = P^/2M* with effective mass M* of the impurity atom [20l[2T]. In the weakly interacting regime 
and sufficiently small momenta the depletion of the host liquid can be studied within the framework of Bogoliubov 
theory. The resulting excitations, polarons, consist of impurities surrounded by a cloud of linear excitations of the 




P/fin 

FIG. 1. Schematic dispersion relation for a mobile impurity in a weakly interacting superfluid. The dispersion defines the lower 
bound of the many-particle excitation spectrum. In the region above the dispersion (shaded gray) there exists a continuum of 
many-body states. The energy is measured in units of the maximum energy of soliton excitations in the superfluid background, 
shown by the dashed line. 
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condensate and were studied in Ref. [32] and later in Refs. [231 [21] ■ However, even for a weakly interacting liquid, 
the induced depletion becomes essentially non-linear at larger momenta and/or strong impurity- liquid coupling. This 
brings a qualitative change of the dispersion relation and we introduce a quasiparticle which we call depleton describing 
the impurity dressed by the non-linear depletion cloud. 

The remarkable property of depleton dispersion is that it is a periodic function of the momentum P with the period 
2Trhn, where n is the density of Id host liquid. To explain this feature it is worth noticing that, being quadratic at 
small momenta, the depleton energy is less than any of the excitations of the host liquid (see Fig. [T]) . Indeed, the 
low-energy excitations of the host liquid have a sound- like nature with the energy cP > P^/2M*, where c is the speed 
of sound. Therefore the dressed impurity excitation provides the cheapest way for the system to accommodate a small 
momentum. This means that the depleton dispersion relation is defined as the lowest possible many-body excitation 
energy of the system with a given momentum. An example of such dispersion is shown in Fig. [T] Above the depleton 
energy there is a continuum of many-body excitations comprised of the moving impurity and a certain number of 
phonons. This spectral edge is characterized by the power-low singularities of zero-temperature correlation functions 
and was discussed in Refs. [^ [35H?7| . 

One may argue now that in the infinite system the ground-state energy with a given momentum is a periodic function 
of the latter with the period 2TThn. Indeed, it is easy to see that the ground-state energy with momentum 2i:hn vanishes 
in the thermodynamic limit. To this end consider a ring of length L where the spectrum of the momentum operator is 
quantized in units of 2-Kh/ L. If the momentum of each particle is boosted by one quantized unit, the total momentum 
of the system is 2TThn, while the total energy vanishes thanks to the total mass diverging in the thermodynamic limit. 
We thus conclude that the ground-state for a given momentum P which is the dispersion relation E{P) of the dressed 
impurity is a periodic function of momentum. Explicit examples are provided by exactly solvable models [25j and 
weakly interacting Bose-liquid corresponding to Fig. [T] and considered in details in Appendix [C[ 

The physics behind the periodic dispersion relation is in the transfer of momentum from the accelerated impurity 
to the supercurrent in the background liquid: similarly to the density depletion, the moving impurity creates the 
sharp phase drop $ across it. To satisfy the periodic boundary conditions the rest of the liquid must sustain the 
phase gradient resulting in the supercurrent which carries momentum Hn^. While the supercurrent is absorbing 
momentum, it does not contribute to energy. Indeed, as it was already mentioned, in the thermodynamic limit the 
bulk of the liquid is infinitely heavy and thus can accommodate any momentum at no energy cost. As a result, the 
energy and momentum of the depleton core, being periodic functions of the phase drop $ with the period 27r, oscillate 
as functions of the total momentum with the period 2nhn, while the rest of the momentum goes into the supercurrent. 
A similar periodic dispersion relation would arise if the host liquid were considered as a rigid crystal with the lattice 
constant a = 1/n. Then the momentum interval between P = —tnm and P — hirn is nothing but the Brillouin zone 
of such a crystal, while the impurity dispersion, discussed above, is its lowest Bloch band. 

Although the above considerations completely disregard the absence of the true long-range order in the Id liquid 
(either superfluid or crystalline) and thus should be taken with care, the periodicity of the dispersion suggests the 
possibility of Bloch oscillations of the impurity atom subject to an external force [2]. Indeed, if an external force F 
is applied to the impurity atom {e.g. electric field is acting on an ion) the total momentum of the system changes 
linearly with time, P = Ft. For infinitesimal force this leads to adiabatic change of the energy E{P) and the velocity 
V = dE/dP of the depleton which become periodic functions of time with the period 2Trhn/F. It is quite remarkable 
that in such a process, on average, the impurity does not accelerate; moreover, it does not even move. Instead, it 
channels the momentum into the collective motion, i.e., the supercurrent, of the liquid and in the process oscillates 
around a fixed location. As a result, no energy is transferred, on average, to the system from the external potential. 

This spectacular phenomenon is present at zero temperature and under an infinitesimal external force. Both finite 
temperature and a finite force complicate the picture in a substantial way. The aim of this paper is to clarify the 
influence of these two factors on the observability of the Bloch oscillations. In brief, our conclusions are as follows: 
at a sufficiently low temperature there is a parametrically wide range of external forces F„^in{T) < F < i^maxj where 
the Bloch oscillations are observable. Contrary to the adiabatic picture, they are accompanied by a drift and exhibit 
certain amplitude and period renormalization. 

The drift manifests itself in the appearance of an average velocity Vd , superimposed on top of the periodic Bloch 
oscillations. It is a linear function of force at small forces, allowing to define the impurity mobility a, as Vd = <jF . As a 
result, the total energy of the system increases (i.e. the system is heated) with the rate FVd = crF^ . This energy goes 
to the emitted long wavelength phonons, which run away from the impurity with the sound velocity c. The maximal 
force can then be estimated as F^ax ~ c/cr. At a larger force the drift velocity exceeds c, leading to Cherenkov 
radiation of phonons. The phonons take a substantial part of the momentum and thus ruin the Bloch oscillation 
mechanism, discussed above. The impurity motion is then either incoherent drift, or an unlimited acceleration, 
depending on the parameters. 

We derive the exact analytic expression for the drift mobility a expressed in terms of the equilibrium dispersion 
relation of the impurity. It is worth noticing that the mobility a is not the linear response property, despite the linear 
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relation between the drift velocity and the external force. Indeed, the dynamics of depleton in this regime is the drift 
superimposed with the essentially non-linear pattern of Bloch oscillations in which the particle explores the entire 
range of impurity momenta and energy. It is therefore not immediately obvious that the mobility a may be expressed 
in terms of the equilibrium properties. Nevertheless we prove that such a relation does exist. 

The true linear response is associated with the thermally induced Landau-Khalatnikov friction force Ff^{V) cx —T'^V, 
which arises due to the Raman scattering of phonons discussed above. It provides the lower limitation on the externally 
applied force -F'inin(T') = ^fr(V'c), where Vc is the maximal equilibrium velocity given by the maximum slope of the 
depleton dispersion, Fig. [l] Indeed, at smaller external forces the velocity saturates and therefore the depleton 
dynamics is confined to the small momenta and Bloch oscillations do not occur. 

The friction coefficient in the small momentum regime was discussed in Ref. |14] and found to be vanishing in exactly 
solvable cases. The reason behind it is the presence of infinite number of the conservation laws, which prevent a non- 
equilibrium state from thermalization. Here we extend those calculations to the case where the impurity explores the 
entire range of momenta and derive an exact result for the full momentum dependence of the friction force. Quite 
naturally, it vanishes in exactly solvable cases too. 

The outline of the paper is as follows: in Section [ll] we give a qualitative description of the depleton quasiparticles 
and the mechanism of coupling to linear sound excitations, or phonons, of the background liquid. Section III is devoted 
to the formal derivation of the depleton Lagrangian based on superfluid thermodynamics. We discuss coupling of the 
depleton with the phonon subsystem and derive a set of stochastic equations of motion for the impurity dynamics 



in Section IV These equation are then used to study radiation losses and derive expression of depleton mobility in 
Section |V] and thermal friction in Section IVII The main results are summarized in Section IVIII Technical details are 
delegated to Appendices. 



II. QUALITATIVE ANALYSIS 

The key to understanding the impurity dynamics is in its interactions with the phonons of the host liquid. This 
problem is rather non-trivial even if the impurity is weakly coupled to the liquid. Indeed, no matter how weak 
the interactions are, the impurity develops local depletion, which become appreciable when the impurity momentum 
approaches im (we set h = 1 throughout the rest of the paper). To visualize this process it is useful to assume a 
semiclassical picture of the background, valid for weakly interacting Bose liquid. In this regime, to accommodate the 
total momentum 7rn the depletion cloud takes the form of the dark soliton [5S]. The dark soliton is an essentially 
non- linear mesoscopic object, which includes a large number of particles and a complete depletion of the liquid density. 
It is exactly the soliton formation which is responsible for channeling momentum into the supercurrent and thus for 
the Bloch oscillations. It is also the soliton which determines the interactions of the impurity with the dynamically 
induced phonons. Therefore the non-equilibrium dynamics of the quantum impurity cannot be separated from the 
dynamics of the essentially non-linear soliton-like depletion cloud. 

What makes the problem analytically tractable is the scale separation between the spatial extent of the local soliton- 
like cloud and the characteristic phonon wavelength. The former is given by the healing length ^ = 1/mc. The latter 
appears to be much longer than ^, if the temperature is sufficiently low and the external force is not too strong. One 
can thus separate the near-field mesoscopic region, which contains the quantum impurity and its depletion cloud, from 
the far-field region, supporting the radiation emitted by the impurity. Since the depletion is restored exponentially 
at the healing length away from the impurity, the precise position of the boundary between the near-field and the 
far-field regions is of no importance. 

Because of the wide difference in their spatial scales, the impurity together with its entire non-linear depletion 
cloud represents a dynamic point-like scatterer for the long wavelength phonons. From the viewpoint of such phonons 
any point scatterer may be entirely described by two phase shifts. These two phase shifts are the discontinuities of 
the phonon displacement and momentum fields across the scatterer. They may be expressed through the number 
of depleted particles N and the phase drop $ across the depleton quasiparticle. Therefore out of many degrees of 
freedom of the near-field region only TV and $ interact with the phononic sub-system. 

What remains is to describe the dynamics of the local depletion cloud with certain fixed values of N and <&. Solution 
of this latter problem is facilitated by the fact that the characteristic equilibration rate of the cloud, estimated as 
1/r = c/^, is much faster than the relevant phononic frequencies. As a result, the cloud may be treated as being in the 
state of local equilibrium, conditioned to certain values of the slow collective coordinates N and $. The fact that there 
are two such slow variables is due to the presence of the two conservation laws: number of particles and momentum. 
The fast internal equilibration of the near-field region is therefore conditioned by the instantaneous values of the two 
conserved quantities. 

The slow change (compared to the fast time scale of r) in the number of depleted particles N and the phase drop 
are due to the fact that the local chemical potential and the local current j at the position of the impurity are 
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both affected by the state of the global phononic sub-system. As a result, one may express the Lagrangian of the 
near-field region as a function of N and through the equilibrium thermodynamic potential which is a function of 
fi and j. This latter function may be independently measured or analytically evaluated in certain limiting cases and 
for exactly solvable models. 

These considerations allow one to separate the local, non-linear but equilibrium problem, from the global, non- 
equilibrium but linear one. The latter statement implies that the host liquid sufficiently far away from the impurity 
may be treated as the linear i.e. Luttinger liquid |29l 130] . This is certainly an approximation which disregards the 
possibility of the moving impurity to emit non-linear excitations, such as grey solitons or shock waves. The train of 
solitons emitted by the impurity moving with a constant supercritical velocity was indeed observed in simulations of 
Ref. [TJ . The kinematics of this process suggests that it is only possible if the drift velocity is close to the speed of 
sound c. We therefore assume that as long as F < Fmax, one may disregard solitons emission and treat the liquid away 
from the impurity as the linear one. This is essentially the same criterion, which allows us to separate the depletion 
cloud from the long wavelength phonons. 

Adopting these approximations, one is able to integrate out the phononic degrees of freedom, characterizing the 
liquid away from the impurity. It reduces the problem to the dynamics of the impurity described by its coordinate 
X{t) and momentum P{t) along with the dynamics of its near-field depletion cloud fully described by the two collective 
coordinates i.e. number of depleted particles N{t) and the phase drop We derive an effective action written in 

terms of such an extended set of degrees of freedom. Such an action leads to the coupled system of quantum Langevin 
equations governing the dynamics of the depleton. 

Away from equilibrium, for F > i^min, the equations of motion yield the pattern of Bloch oscillations. The 
deterministic part of these equations provides with the information about drift velocity, amplitude and shape of the 
velocity oscillations, as well as their period. The stochastic part results in a certain dephasing of the oscillations. 
It is interesting to notice that the stochastic part is manifestly different from the equilibrium noise, prescribed by 
the fluctuation-dissipation theorem. As a consequence, the exactly integrable models loose their special status and 
their non-equilibrium dynamics appears to be not qualitatively different from the dynamics of generic non-integrable 
models. 



III. LAGRANGIAN OF THE MOBILE IMPURITY 



Let us first consider the background liquid in the absence of impurity employing hydrodynamical description 
proposed by Popov |29| . Its Lagrangian is expressed in terms of the slowly varying chemical potential /i and density 
n as an integral of the local thermodynamic pressure 



Lo(/i,n)= dxpQ{iJ.,n) — I dx fin — eQ{n) . (1) 



Here eo(n) is the energy density of the liquid. In the thermodynamic equilibrium the density is a function of the 
chemical potential, given by a solution of the following equation: jj, = ii{n) = deo/dn, which is a result of the 
minimization of this functional with respect to n. This way one defines the grandcanonical thermodynamic potential 
of the host liquid as 

no{fi)^-Lo{ti,n{fi)). (2) 

For a uniform system the Lagrangian and the corresponding thermodynamic potential are both proportional to the 
length of the system. 

Consider now an impurity of mass M, having a coordinate X and moving through the liquid with velocity V = X, 
as measured in the laboratory reference frame. It is convenient to choose the reference frame where the impurity 
is at rest and the liquid flows with the velocity —V, as shown in Fig. [2j In this co-moving frame the impurity 
experiences the supercurrent j' and the chemical potential fi' . Hereafter primes denote physical quantities defined 
in the co-moving frame to distinguish them from the corresponding quantities in the laboratory frame. We employ 
Galilean transformation into the moving frame which gives 

j' = -nV, ^i' = ^JL + mV'^/2. (3) 

Together with the Galilean transformation of the energy density e'^ — + mI/^/2, the transformation ([3]) combined 
with Eqs. ([T) and ([2| show the invariance of the background grandcanonical potential i^Q(M') = As expected 

from the Ga ilean invariance, the latter is independent of the velocity V or the supercurrent j' . 

We introduce now the impurity into the fiowing liquid maintaining both j' and /i' fixed and let it equilibrate. 
Its motion distorts the host liquid density and velocity fields, forming the depletion cloud moving along with the 
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impurity. The grandcanonical potential increases by an amount = E'^ — ^'N^, where E'^ and are the 

corresponding changes in energy and number of particles. Using the Galilean invariance one can relate the energy 
E'^ = E^ — P^V + to the energy E^ and momentum induced by the moving impurity in the laboratory 

frame. This fact and relations ([3| allows one to identify, in the spirit of Popov's approach, the Lagrangian of the 
depletion cloud with the negative change of the grandcanonical potential 

Ld(F, n) = P^V -Ea + /iA^d = -E'^ + //^d - m') ■ (4) 

This relation is quite remarkable as the left hand side describes dynamics of the polarization cloud moving with the 
velocity V ^ while its right hand side is the thermodynamic quantity. The link between them comes from the Galilean 
transformation, Eqs. ([3|. The relation between the Lagrangian and the grandcanonical potential, Eq. Q, can be 
viewed as a generalization of the Popov relation, Eq. ([2|, to the case of mobile impurities. 

Two remarks are in order. First, in assuming the thermodynamic equilibrium at nonzero supercurrent j' flowing 
through the impurity we rely on the superfluidity. Second, we note that the increase in energy E'^^ momentum Pj, 
number of particles and the grandcanonical potential due to the presence of one single impurity are finite size 
corrections to the corresponding extensive quantities. 



A. Collective degrees of freedom of the depleton 

Variations of the so far fixed parameters j' and /i' of the background liquid induce changes in the thermodynamic 
potential of the depletion cloud. It can be written with the help of the corresponding response functions as 

df^^ = $ dj' + iVd^' , $ = c»jvO', N = d^,,Vl'. (5) 

The response to the variation of the chemical potential N = — A^d is identified with the number of particles expelled 
from the liquid by the impurity (hence the minus sign). The response $ to the change of the supercurrent j' is 
the superfluid phase and has no analogy in classical thermodynamics. In the state of the global thermodynamic 
equilibrium both <I> and N are rigidly locked to j' and fi' and, consequently, to V and n. This is denoted by writing 

= 3'o(^, n) and N = Nq{V, n). These functions can be obtained from the derivatives of the Lagrangian defined in 
Eq. Q as described in the next subsection. 

In the nonequilibrium situations, where the supercurrent and chemical potential fluctuate it is convenient to treat 
$ and N as independent variables. We perform the standard Legendre transformation to a new thermodynamic 
potential, 

i/d($,iV) = r!^-/$-Ai'iV, diJd = -/d$ - //'diV. (6) 

The independent variables $ and N describe the state of the depleted liquid in the immediate vicinity of the impurity 
which may or may not be in equilibrium with globally imposed j' and )j! . In equilibrium situation Ha{^, N) does not 
contain any additional information with respect to the thermodynamic potential , fi'), which is in turn related 
to the depletion cloud Lagrangian Ld{V,n) by Eq. Q. The aim of introducing Hd{^,N) is to allow for interactions 
of the depletion cloud with the long wavelength phonons. As we shall see in Section |IV| the latter may change the 
number of particles and the momentum of the depletion cloud, forcing it to equilibrate to some new values of and 
A^. Before turning to phonons, it is instructive to rewrite down the impurity Lagrangian Q by substituting into it 
Eq. (|6| and considering <I> and A^ as independent vaiables. 



Ld = - MV^ - /$ - fi'N - i7d($, N) . 



(7) 
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Expressing j', ^' through V , n using Eq. (|3|, we finally obtain the Lagrangian of the depleton 

L{V, $, TV) = ^ {M - mN)V^ + nV^ - - i?d($, N) , 



(8) 



where = fj,{n) is the equilibrium chemical potential of the host liquid in the laboratory frame and we have added 
the bare kinetic energy proportional to the impurity mass M. The momentum of the depleton is obtained by the 
standard procedure 



P 



(M - mN)V + n$ . 



(9) 



The last term describes the supercurrent momentum n$ stored in the background. The first term is proportional to 
the reduced mass M — mN, expressing the fact that the mass mN is removed from the local vicinity of the moving 
impurity. This quantity should not be confused with the effective mass M* given by the curvature the equilibrium 
dispersion Eq. (15) or obtained from the equilibrium dynamics in Eq. (19 1. 

We can use Eq. ^ to express velocity of the depleton as a function of its momentum, 



X = V ^ V{P, $, N) 



P -n^ 
M - mN ■ 



Combining Eq. ( 10 ) and the Lagrangian ([8|) leads to the Hamiltonian 

1 CP — n^)'^ 

This Hamiltonian generates the following equations of motion 

P = 

= -dg.H ^ nV - dg.Hd 
= -OnH = -/i - my 2/2 - 



(10) 



(11) 



(12) 
(13) 
(14) 



in addition to Eq. (10). The first equation is the momentum conservation expected in the absence of external forces 
and for homogeneous background. Two other equations are in fact static constraints. This is a manifestation of the 
already mentioned fact that without phonons 4> and N are rigidly locked constants and do not have an independent 
dynamics. 



B. Dynamics of depleton in the absence of phonons and equiHbrium values of collective variables 



As we shall see in Section |IV[ the coordinates canonically conjugated to N are phononic displacement and 
phase at the location of the impurity. In the absence of these degrees of freedom the only consistent solution of 



Eqs. (13), (14) is static relations $ = ^o{P,n), N = No{P,n). Substituting them back into the Hamiltonian (11) 



leads to the equilibrium dispersion relation of the dressed impurity (depleton) , 

H{P, $o(P, n),NoiP, n)) = EiP, n) . (15) 

The corresponding "equilibrium Lagrangian" can be obtained as L(V, n) = PV — E hy expressing the momentum 
with the help oi V — dE/dP. In most situations it is rather these quantities and not the "internal energy" H^{^,N) 
represent the physical input about the dynamics of the dressed impurity. They can be obtained from solving the 
equilibrium problem for the impurity moving with the constant momentum P or velocity V, through the liquid with 
the asymptotic density n. Below we show explicitly how collective variables $ and N and the corresponding energy 
Hd can be obtained from the knowledge of E{P, n) or L(V, n). 

We start with the situation when velocity y is a control parameter. In this case finding Hd{^,N) amounts to 
perform the Legendre transformation Eq. ^ by exploiting the definition (|4| of the Lagrangian = L—MV'^/2 of the 
depletion cloud in terms of the thermodynamic potential ^l'^ and using the pair (V, n) instead of the thermodynamical 
variables To achieve this goal we use Eqs. ^ to relate the corresponding partial derivatives by the linear 

transform 



d„J ^ \-V mc^/n) \ d^, 



(16) 
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Here we have used the relation dfi/dn = mc^/n between the compressibility and the sound velocity c. Using the 
definitions Eqs. ^ and Eqs. (16) we can express the derivatives of L{V,n) in terms of collective variables 



dL_ 
dV 



- MV = P - MV = n$ - mVN 



dL 
dn 



mc 



N. 



(17) 
(18) 



Solving these equations yield equilibrium values, ^o{V,n) and No{V,n). Eq. ( [T7| can be otherwise obtained by 
simply substituting ^Q{V,n) and NQ{V,n) into the definition of the momentum, Eq. (|9|. This is a consequence of 
the equations of motion (Il3|), (14). The quantities involving second derivatives of the Lagrangian do not enjoy this 



property. The most obvious case is the effective mass 



M* 



dP 

dV 



M — mN — mV 



dNo 9$o 



dV 



dV 



(19) 



which differs from the expression M — mN obtained by taking partial derivative of Eq. ([£]) with respect to velocity. 

The equilibrium relations $o(V^, n) and iVo(y,n) can be inverted to find velocity Vb($,A^), and density no{^^N) 
for given values of $ and N . Substituting them into Eq. (Is]) gives 



i/d($,A^)--i("t^o,no) + 



MVg 



N . 



(20) 



Conversely, we can use the momentum P as a control parameter. Again, using equations of motion Eqs (13), (14) 
one is able to show the equivalence of the derivatives dH{P,^,N)/dn — dE{P,n)/dn and dH{P,^,NjJdP = 
dE{P,n)/dP. The latter defines the velocity Vo{P,n). Using these facts and differentiating explicitly Eq. (11) we 
have the following system of equations 



Vo{P,n) 



M - mN ' 



dE{P, n) 
dn 



-iV-Uo(F,n)$, 



(21) 



which are equivalent to Eqs. (17l, (18) by virtue of the fact that {dL/dn)v = —{dE/dn)p. Solving Eqs. (21) 
yield equilibrium values ^q{P, n), Nq (P, n) as functions of P and n. Next, we invert these relations to obtain 
Po{^, N),nQ{^, N) as functions of $,iV. Using Eq. (15) and Eq. (11) we obtain the expression for the core energy 



Hd{^,N) ^ E{Po,no) - ^iino)N - 



1 (Pq ~ no$)- 

2 M-mN 



(22) 



in terms of the dispersion E(P, n) and its derivatives. 

To illustrate this procedure we use two cases, where the energy Hd{^,N) possesses a simple analytic form. One is 
the grey soliton in weakly interacting Bose-Einstein condensate describing a massless, M = 0, impurity propagating in 
a weakly interacting Bose liquid with the coupling constant g. The standard results [^51 [31] for the soliton dynamics 
are provided in Appendix [B] and lead to 



(23) 



Another example is provided by a strongly interacting impurity |211 132j . In this case the number of expelled particles 
TV is almost independent on the state of the impurity and may be considered as a non-dynamic constant. The 
remaining dependence of energy on the superfiuid phase $ has a standard Josephson form 



i7d($) = -Pjcos$. 



(24) 



The Josephson energy Ej — hnVc is expressed through the corresponding critical velocity 14, which in this case is 
much smaller than the sound velocity c. 



Expressions ( 20 ) or ( 22 ) provide the core energy of the locally equilibrium depletion cloud as a function of its slow 



variables <i> and N. This procedure emphasizes the fact that introduction of $ and iV does not rely on the semiclassical 
interpretation of the condensate wavefunction. In fact, they may be defined even away from the semiclassical weakly 
interacting regime, where the phase of the condensate as well as its depletion are not well defined. 
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IV. COUPLING TO PHONONS 



An external force F acting on the impurity drives the system away from equihbrium, making the impurity radiate 
energy and momentum. For a sufficiently weak force such a radiation takes the form of long wavelength phonons, 
i.e. small deviations of density p{x,t) and velocity u(x,t) fields from their equilibrium values, see Fig.[3| Below we 
show how coupling to phonons can be formulated in terms of the collective variables ^,N. It turns out that this 
procedure is based solely on the principles of gauge and Galilean invariance, i.e. conservation of number of particles 
and momentum, and leads to universal results. 



A. Hydrodynamic description of phononic bath 



We start by considering the Lagrangian, governing dynamics of the phonon fields in the bulk of the liquid. To 
this end it is convenient to parameterize them by introducing the superfluid phase ip{x,t) and the displacement field 
"dix, t) such that u = dxf/m and p ~ OxiS/tt. The dynamics of these variables can be described following the method 
of Popov [29] by considering slow change of the density n — > n + p{x,t) and, independently, the change of the chemical 
potential, 



fj.^ p- ip{x,t) - 



mu^{x, t) 



Substituting them into Eq. ([T]) yields the Lagrangian of phonons. 



dx 



Paip{x, t),nix, t)) - po{p, n) 



dx 



-pip- 



(^eo(n + p)~ eo(n) - pp'j 



(25) 



(26) 



For nonzero phononic fields, the impurity is subject to the modified local supercurrcnt and chemical potential in 
the co-moving reference frame 



p =p-(p 



m{V 



J 



2 

P){V- 



u) 



2 

-nV 



p{n) 



{if + Vdx(p) ; 



(27) 
(28) 



where the phonon variables are taken at the instantaneous spatial position, X{t) of the impurity. Equation (27) 
follows from Eq. (25) and the fact that in the presence of the background flow u = dx^/m the velocity of impurity 
with respect to the liquid is changed to V — u. To derive expression (28) for the modified supercurrent we have used 
the continuity equation in the form "iJ/Tr — — (n + p)u. This relation is an exact statement, which follows from the 
gauge invariance and is valid for any configuration of the fields. 

Substituting the modified supercurrent and chemical potential, Eqs. (27), (28), into Eq. ([t]) and subtracting the 



corresponding equilibrium values, results in the following universal form of the interaction Lagrangian, 



dt' 



(29) 



It is full time derivative d/dt = dt+Xdx — dt + Vdx which enters the interaction term, as follows from Eqs. (27), (28 1 



The interaction Lagrangian, Eq. (29) provides dynamics of the collective coordinates N and $. It shows that the 



corresponding canonical momenta are the phonon degrees of freedom at the location of the impurity, i.e ip{X,t) and 
'd{X, t) correspondingly. Through the gradient terms in Eq. ( 26 ) these two local variables are connected to the phonon 



fields elsewhere and it is the dynamical properties of these phonons which determine the behavior of the impurity. For 
example if the spectrum of phonons is discrete, one expects coherent oscillations of few modes. In the infinite system 
the continuous spectrum of the background modes leads to dissipation similar to that of Caldeira-Leggett model [33] . 




FIG. 3. Impurity propagating in local environment. 
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B. Linear phonons and transformation to chiral fields 



Away from this interaction region the excitations of the hquid may be considered as Hnear ones. Therefore one 
may keep only the quadratic terms in the phononic Lagrangian Eq. ( 26 ) . For the kinetic energy in Eq. ( 26 ) one thus 
retains the leading term mnu^/2, while the potential energy is expanded as eo{n + p) — eo{n) — fip {dJIJdn) p'^ / 2 = 
(mc^ /2n)p^ , using thermodynamic relation between the compressibility and the sound velocity c. As a result, one 
obtains the quadratic Luttinger liquid Lagrangian 



Lph ^ - I dx 

TT 



(30) 



Besides the sound velocity c the Lagrangian in Eq. ( 30 ) is characterized by the dimensionless Luttinger parameter 
K = im/mc, which depends on the degree of correlations in the host liquid [3^. For a liquid of weakly repulsive bosons 
the Luttinger parameter K is large, K ^ 1 but can be reduced down to the limiting value K — Ihy increasing the 
repulsive interactions or decreasing the density of the particles |34| . The combination irc/K = mc^ /n entering Eq. (30) 
contains information about interactions between the particles in the liquid, while the combination cK/tt — ck — n/m 
is independent of the interactions as a consequence of the Galilean invariance. 

To gain an additional insight into the physics of the depleton-phonon interaction, the phononic fields can be 
decomposed into a doublet of right- and left-moving chiral components with the help of the linear transformation 



?9/7r 



X- 



T 



1 

\/2 ^7^/ 



(31) 



where k = K/t: = n/mc. In terms of the chiral fields the Lagrangian, Eq. (30) splits into a sum of two independent 
contributions. 



^phM = ^ j dx [x+{d.^dt + cdl)x+ + X-{-dxdt + cd^X- 
The matrix of inverse phonon propagator is defined by its Fourier representation. 



^ldxx^D~\. 



S{x)S{t)D- 



dq duj 
27r 27r 



q{uj - cq) 

— g(i^ + cq) 



(32) 



(33) 



The equation of motion following from the Lagrangian (32) dictate a simple coordinate and time dependence, 



X±(a;,t) — x±i^ T ct). Using this fact one can show that for imiformly moving reference point X — Vt one has 

d 



X±{X,t) = {VTc)dxX±{X,t). 



(34) 



This property is in fact a statement about correlation functions of the fields x± calculated with the Gaussian action, 
Eq. (32) which enforces classical equations of motion; the path integration is performed over arbitrary configurations 



of the fields. The interaction term, Eq. (29) may be conveniently rewritten by introducing chiral collective variables 



A = 



(35) 



In the static limit the quantities K± are proportional, up to factor of to the chiral phase shifts 6± introduced in 
Refs. |26l [35] . In presence of external force F they acquire the dynamics, which is governed by the total Lagrangian 
of the depleton interacting with the phonons 



Ltot = PX- H{P, A) + UiX)- At (t) x{X, t) + iph [x] 



(36) 



where the impurity Hamiltonian H{P,A) is obtained from Eq. (|ll[) by using the linear relation, Eq. (35) between 



the chiral phase shifts A and collective variables $ and N, and iph is given by Eq. (32). Here U{X) is an external 
potential acting on the impurity atom only and F = —dU jdX is the external force. 



C. Integrating out the phonons 



We have now all necessary ingredients for describing dynamics of impurity coupled through the interaction term, 
Eq. ( 29 ) to the phononic bath. The presence of the impurity is felt by phonons through time-dependent boundary 
conditions at a; = X(t) parameterized by the collective variables X[t) and ^(t), or equivalently by chiral phases A±(t). 
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Here we simplify even further our description by solving phononic linear equations of motion for any variation of these 
collective variables and substituting the obtained solution back into the action. This procedure leads to dynamics of 
the impurity expressed in terms of collective variables only and is equivalent to exact integration of Gaussian phononic 
action. 

To this end we employ the Keldysh formalism [M1I57| and extend the dynamical variables X{t), P{t), A±{t) as well 
as phononic fields x±{^it) to forward and backward parts of the closed time contour t ^ t±. Performing Keldysh 
rotation, we write them as X{t±) ~ Xci±Xq/2, h.{t±) — Aci±Aq/2, and x(i±) = XciiXq/2 with the help of symmetric 
("classical", cl) and antisymmetric ("quantum", q) combinations. The coupling term Eq. (29) then becomes 



Li, 



-A^l(t)Xq(Xcl,t) - At(i)xel(Xci,<) - A^l(t)X,(t)9,Xcl(Xcl,i) 



(37) 



up to terms linear in Xq. The advantage of chiral fields introduced in Eq. (31) is that one can use the property (34) 
together with classical trajectory X^x = Vt to simplify the interactions, Eq. (37 1 as 



-A! 



At + XqAtjV-l) diXcl 



(38) 



where we have introduced the matrix 



1 

c-V 

- 



c+V 



(39) 



The interaction term, Eq. ( 38 ) is linear in phononic fields so that a Gaussian integration with quadratic action, 
Eq. (36) can be performed by standard methods as explained in Appendix [A] It leads to quadratic, though a nonlocal 
in time effective action for the collective variables. 



5cfr = -2 /d<At;(i) 



Aq(i)+V-iA,i(t)Xq(i) 



Am' 



At (i) + Xq(i)Atj(i)V-ll dtf{t - t') [Aq(t') + V-lAcl(i')^q(i') 



(40) 



where, assuming thermal equilibrium of the phononic subsystem, the matrix F(<) is related by inverse Fourier transform 
to the matrix 



, , /coth^ 
F(^) = 







coth 



(41) 



2T_ 



of thermal distribution of the chiral bosons with the temperatures T± = T{1 ^ V/c) modified by the corresponding 
Doppler shifts. 



One should supplement the action Eq. ( 40 ) with the Keldysh analogue of the action corresponding to the depleton 
Hamiltonian, Eq. (Ill, 



S^ldt {X,i-dpH]P^-[P,i-dxU)X^-VAH -A^ , 



(42) 



where we kept only terms linear in the quantum components and H — H{Pci, A^i). Notice that quadratic terms in 
quantum fields are absent in Eq. ( 42 ) while cubic and higher orders are omitted in the spirit of the semiclassical 
approximation. 



The second line in the effective action Eq. ( 40 ) may be split with the help of the Hubbard-Stratonovich transforma- 



tion, which introduces two real uncorrelated Gaussian noises ^+(/:),^_(t). Their correlation matrix in the frequency 
representation takes the standard Ohmic form (see, e.g. [33]), 



e(c.)et(c.)\=c.F(c.) 



(43) 



The action, Eq. (40) becomes local in time, 

SeS^-lJdt (At, -2^1 



Aq+V-lAelXq) 



(44) 



12 



Now the entire semiclassical action is linear in quantum components and integration over them enforces the delta- 
functions of the equation of motions. While Eq. (10) remains intact, due to the absence of Pq in the effective action 



(|44|), Eqs. ^l2\, and ^T^ are modified by the phonons 



(45) 
(46) 



where we have dropped subscripts for clarity. The obtained equations include additional dissipative terms involving 
time derivatives of the collective variables A. They also include fluctuations coming from the pair of Gaussian noises 



S,±(t) correlated accordingly to Eq. (43). 



DEPLETON DYNAMICS AT ZERO TEMPERATURE. RADIATIVE CORRECTIONS 



Our goal is to discuss non-equilibrium solutions of the equations of motion in the presence of a constant external 
force F. Neglecting for a moment fluctuation terms and using transformation Eq. (35) the equations of motion, 



Eqs. (45), (46) can be rewritten in terms of the collective variables $ and N as follows, 

'kV 



P^^-^(<i>,iv)TV-Tt [l 



= F - 



c? - V"^ 



2c 



V 

2kc' 



TV 

2k 



— nV - 9$-ffd 



(47) 

(48) 
(49) 



The rate of energy radiated by phonons is obtained by taking derivative with respect to time of the total impurity 
energy. 



W = H - FX ^ V{P - F) - • VaF = V{P ~ F) + ^dg.H + NOnH. 
Using Eq. (10) and equations of motion either in the form (45^46) or (47-49) we obtain 



W : 



-A^ [i + vy-^]A = -- 



y2 



V 



2k 



(50) 



(51) 



The dissipation of momentum, Eq. (47) and energy, Eq. (51) is a generahzation of Eqs. (22), (23) in Ref.[38] (up to a 
factor of two), where they were derived in the context of dynamics of grey solitons. 

According to Eqs. (12)-(14), in the absence of the external force F — Q there is a family of stationary solutions 
of the equations of motion, which are characterized by a constant velocity V below some critical velocity Vc- This 
solutions describe the dissipationless motion of the impurity consistent with the superfluidity. Indeed, by neglecting 
the fluctuation terms we effectively put the temperature to be zero, thus making the one-dimensional liquid superfluid. 

On the first glance one can just solve the set of the evolutionary equations (10), (47) -(49) to fully describe the 
impurities dynamics. One needs to be careful, though, because Eqs. (48 ), (49 ) correspond to the motion in the vicinity 



of the maximum of the Hamiltonian H and therefore exhibit runaway instability. The further look at this instability 
shows that its characteristic rate is of the order mc^ ^ ^, which is well outside the frequency range of applicability 
of the theory developed above. In fact this high-frequency instability of $ and N evolution is a direct analog of 
the well-known spurious self-acceleration of charges due to the back reaction of the electromagnetic field [39]. The 
recipe to overcome it is, of course, also well-known: instead of trying to solve equations of motion directly, one should 
perturbatively find how radiation corrections influence the dynami cs |40| . This strategy offers a convenient analytical 
approach to treat the dynamics described by Eqs. (10) and (47)-(49). Below we apply it to study modifications to 



Bloch oscillations which arise due to the phonon radiation 



A. Radiative corrections to Bloch oscillations 



At zero temperature with a sufficiently small applied force one expects the system to adiabatically stay in the 
ground state with total momentum P, while the latter is changed by the external force P = Ft. In this zeroth 
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approximation, the motion is nothing but a tracing of the dispersion relation E{P, n) and the phononic subsystem 
gains no share of the work done on the system by F. The velocity is simply 



= VQiFt,n) = dE{P,n)/dP 



P=Ft 



(52) 



Since the dispersion relation displays a periodic behavior, one immediately obtains velocity Bloch oscillations with 
period t'^^ — 2TTn/ F^ amplitude and zero drift. In reality, the slow acceleration of the impurity over the course of 
a Bloch cycle gives rise to a soft radiation of low energy phonons, which serve to renormalize the period, amplitude 
and drift from the zeroth order approximations. 

To study the corrections to the dcplcton trajectory, let us assume it exhibits a steady-state motion such that 
T/(i+TB) ^V{t), N{t+TYi) ^ N{t) and<I>(t+rB) = $(i) + 27r. Then it follows from Eq. ([9| that P(t+TB) = P{t) + 2Tm. 
Here tb is the, a priori unknown, true period of the motion, not to be confused with the zeroth approximation Tg*^^ 
To find it we integrate Eq. ( 47 ) over a single Bloch cycle 

1 



27rn 



dtp = Fn 



B 



dt- 



1 



,2 _ y2 



(53) 



where -Frad is the average radiative friction force exerted on the impurity over a single Bloch cycle. Since the radiative 
frictional force tends to reduce the applied force, Eq. (53) indicates that the true period of oscillation is larger than 
the zeroth approximation 



(0) 



The work of the external force per unit time is given by 

FV = E[P,n) 



W. 



(54) 



The first term in the r.h.s. of this equation is the reversible change in energy of the impurity, while the second term, 
owing to Eq. (51), is the rate of energy channeled into the phonon system. We average Eq. (54) over a single Bloch 



cycle, noticing that {E) = due to the periodicity of the dispersion relation. The remaining term corresponds to the 
power radiated into phononic bath and leads to the drift velocity: 



-{W)/F 



(55) 



Assuming the energy pumped into the phonon system per Bloch cycle to be small we use the bare trajectories. 
Vo{P, n), ^oiPi ^oiP, n). Using the fact that dt = dP/F one can show that 

Fd = crF. (56) 



Here a is the T 



mobility of the impurity, given by the average over the Brillouin zone 

2 



1 

27rrt 



dP 



1/2 
•^0 



dP 



c 



dP 



dP 



1 
2k 



dP 



(57) 



It was mentioned in Sec. IIIB that the equilibrium functions <I>o(^', and NQ{P,n) can be obtained directly from 
partial derivatives of the equilibrium dispersion relation E(P,n). Since Vq = dE/dP, the mobility may be expressed 
entirely in terms of E{P, n) and the Luttinger parameter K = ttk. 

The fact that the mobility can be expressed through equilibrium properties is reminiscent to the Kubo linear 
response formulation. It is crucial to mention that they are not a linear response property in the Kubo sense. The 
Kubo linear response takes place at finite temperature, see Section VI At T = the liquid is superfluid and the 
impurity undergoes Bloch oscillations with the amplitude Vc at arbitrarily small external force F. It means that 
the response is essentially non-linear. The mobility a describes the average (over one period) shift of the oscillation 
center due to the energy radiated in the course of such non-linear oscillations. The fact that it may be fully expressed 
through the equilibrium properties is rather remarkable on its own right. 

The result in Eq. ( 56 ) holds for sufficiently weak external perturbation F 



< F„ 



comparing the corresponding drift velocity with the velocity of sound. From Eq. (56 1 we have 



One can estimate -Finax by 



-Flnax = c/ a. 



(58) 



As the force increases the separation of length and energy scales used to define the depleton dynamics cannot be 
justified. In other words for a strongly perturbed system the equilibrium dispersion relation E{P,n) ceases to be a 
meaningful concept. 

Below we use a model of an impurity coupled via delta-function interaction with strength G to the background 
particles to illustrate the dynamical properties of a depleton discussed above. We discuss two regimes: the strong 
coupling regime, where the interactions in the liquid can be of arbitrary strength and the weak coupling regime where 
the dynamics of the depletion cloud is governed by Gross-Pitaevskii equation. In the latter case one is restricted to a 
weakly interacting bosonic background. 
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FIG. 4. Velocity as a function of time for various forces listed in the legend {F in units of Fmax ~ 2nmc^). Here k = 20, 
G/c — 20 and M = 40m. The dashed lines correspond to the drift velocity plotted in Fig. |5] One notices that as F increases, 
the drift velocity and frequency of oscillations increases, while the velocity amplitude (as measured from Vu) decreases. 



1. Strong coupling regime G/c ^ 1 



The impurity expels a large number N ^ 1 of particles from its vicinity thus one can neglect the dynamics of N 
and use the Josephson form, Eq. (24 1 for the remaining dependence of energy on superfluid phase. Due to the fact 
that velocity is bound by the critical val ue Vc the impurity is slow V < Vc <^ c. Indeed, the calculation for weakly 



interacting bosonic background, Eq. (C9) (see also Ref.[41) ) gives Vc = c /2G <^ c. Therefore, if the impurity is not 
too heavy, M ^ 
contribution P 



le second term in the momentum, Eq. ^ and we are left with the superfluid 



m, we can neglect t : 

= n<E> only. Using this fact in Eq. (|57|) gives the universal result for the mobility 



K 



K _ 1 
h-n? 2nmc ' 



(59) 



where we restored h = h/2'n. Note that in the strong coupling limit, the mobility is independent of the impurity 
parameters and only depends on the parameters of the host liquid, namely the Luttinger parameter K and the 
asymptotic density n. This result has been obtained by Castro-Neto and Fisher, [TJ] by using the linear response 
approach. 

For the case of impenetrable bosons or free fermions corresponding to K — \ one can use the analogy from the 
electronic transport. Suppose the background is made of non- interacting fermions each carrying electric charge e. In 
the frame co-moving with the impurity a current I = enVi) flows through the wire, whose quantum resistance is R. 
The latter is given by the Landauer formula, which for spinless case reads as R — h/e^ . The ohmic power transferred 
to the system, /^i?, must be supplied by the external force, FVd = hn^V^ giving the result in Eq. ( |59| for K = 1. 
For K > 1 one recalls that R = h/Ke^ |16j leading again to Eq. (59). Notice that discussion of Ref~[42) claiming 



interaction-independent mobility is not applicable here, since we always assume that the system length is much larger 
than the characteristic wavelength of phonons. 

Scaling with K implies that the mobility is higher for weaker interactions and diverges in the free-boson limit. To 
understand this result, notice that the moving impenetrable impurity experiences nVu collisions per unit time. A 



given collision results in the momentum transfer toVd to the boson of the gas. The balance of forces F = nmV^ then 



implies Vd = F/\/ nmF leading to u — 1/^/nmF — )■ oo as — )■ 0, in agreement with non-interacting limit of Eq. (59 1. 

Turning to the period of oscillations, Eq. (53) we use the relations V ~ V-q = nF/2n'^, <I> w F/n to obtain the 
renormalized period of Bloch oscillations, 



2'Kn 
~F~ 



( 



(60) 



The characteristic force entering Eq. ( 60 ) coincides with the upper bound i^max — 2nm<? obtained from Eq. ( 58 ) 
using the result (59) for the mobility. Indeed, the frequency of the motion oj « F/n should be small compared to the 
typical phonon frequency mc^ in order to justify the large scale separation employed in Sec. H 
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FIG. 5. Drift velocity Vd (left panel) and amplitude Vb (right panel) as functions of the applied force F corresponding to the 
same set of parameters are as in Fig. [4] The solid curve on left panel is indistinguishable from the analytic prediction, the first 
of Eqs. (641, while the dashed curve is the small force result Eq. (56 1 with mobility a — K/2hn^. On right panel the solid curve 
is numerics, the dashed curve is Eq. (|D9|, and the dotted curve corresponds to the small F limit of Eq. (65 1. 



To go beyond lowest order in F we devise a numerical approach to the strongly coupled impurity based on the 
Josephson form, Eq. (24) for the energy. Using it in Eq. (48) and recalling Eq. (47) gives, 



P = F ~ - 



1 nV 



2 c2 - T/2 



= V -Vr sin$, 



(61) 



where P is given by Eq. ([9]). The second of Eqs. (61) relates the change in $ with the deviation of the impurity 
velocity V from its equilibrium value VcSin<i>. The velocity V{t) is obtained by solving Eqs. (61) numerically and the 
results depicted in Figs. |4] agree well with the ansatz ( 62 ) . As one can see the precise choice of initial conditions is 
essentially irrelevant to the ensuing discussion, where we focus only on the asymptotic, steady state behavior of the 
functions $(i) and V{t). 

After a sufficiently long time the system reaches the regime of steady Bloch oscillations, where its velocity is given 

by 



V{t) = Vt) + Vb cos{ujBt + 5) . 



(62) 



Here the parameters Vd, ^b, "^b — Stt/tb and S depend on the ratio of the external force F to the critical force Fmax- 
To find the drift velocity and period of Bloch oscillations we substitute the ansatz V = Vd and P = ni> = ojb into 



Eqs. (61 ) to obtain a closed set of equation. 



^ _ 1 kVdujI^ 



2c2 



1/2 ■ 



LOB 

2mc2 



which are solved to give 



c 



_Wb_ 

2mc2 




(63) 



(64) 



Expanding the expression for ojb ~ 27r/TB for small F/F-n-iax one recovers the result (60) obtained by a different 
method. The results for the drift velocity are plotted in Fig. [5^ and are in excellent agreement with numerical 
solution of Eqs. (61 ). 



The Bloch amplitude Vb can also be estimated from Eqs. ( [6l| ). In Appendix [d| it is shown that 



Vb = VA\ 



1 



M - mN 
rriK 



F 



(65) 



Thus, Eq. (65) predicts a decrease of the Bloch amplitude Vb with increasing F (see Fig. [5p). Physically sensible, 
this result implies that as the force increases, the ideal equilibrium tracing of the dispersion relation is, to a degree, 
lost. This is due to the phononic subsystem gaining a proportionately larger share of the work done by the external 
force. Expression (65 1 is compared with results of the numerical solution of Eqs. (61 ) and the results presented in the 
right panel of Fig. 5p show a good agreement. 
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2. Weakly interacting background bosons and weak coupling regime G/c<^ 1 



We deal with the case of the background made of weakly interacting bosons by using the Gross-Pitaevskii equation 
for the background Bose liquid as explained in Appendix [C] By calculating numerically the dispersion relation of the 
impurity, see e.g. Fig. [l] we obtain the functions A^o and $o for all values of the impurity-background coupling G. 
The results of numerical evaluation of Eq. ( 57 1 are shown in Fig. |6] 

Let us turn now to the case of weak coupling G c, where particles can tunnel through the semi-transparent 
impurity. In such a case one expects the mobility to increase as the impurity becomes more transparent. For weak 
coupling, the main contribution to the integral in Eq. (57) comes from the region where the velocity of the impurity is 
maximal, i.e. V ^ Vc- It corresponds to the inflection points of the dispersion, Fig.[l] Beyond this point the dispersion 
follows closely the dispersion of grey solitons, Eq. (B4| and one can use the expressions (B5) to estimate the mobility. 



As the momentum Pc corresponding to maximum velocity V{Pc) = Vc is small we simplify Eqs. (B4|, (B5| and write 



1 /$ 

2 V 2" 



$ - sin $ 



$3 



(66) 



or, equivalently. 



$(P) = (6P/n)i/3; 



V{P) 



{6P/n) 



2/3 



(67) 



Substituting these expressions into Eq. (57) we can estimate the mobility as 

dP/n K 



K j dP (d±\ 
'^^ njp^ l~Vyc^ \dP ) 



K 



i2 G 



1/nmG . 



(68) 



Here we have estimated the momentum cutoff A — Pc/n 
equation in (67) 
M 



G/c using Eq. ( |C7[ ) for the critical velocity and the second 
The coefficient in Eq. (|68| is fixed from the numerics and was found to be very close to one for 



m. Interestingly, in the weak coup 



ing limit a is independent of the correlations within the liquid (provided 
the liquid is weakly interacting, G <C c), diverging in the limit of a completely transparent impurity. Thus, in the 
weak coupling limit, the upper critical force F^ax = nGmc corresponds to energy difference Gn per healing length 
^ ~ 1/mc. 




G/c 



FIG. 6. Log-log plot of the zero temperature mobility ct as a function of the impurity coupling G in a weakly interacting 
superflui d (s olid li ne) . The dashed lines correspond to the asymptotic values of a in the strong and weak coupling limits given 
by Eqs. (59 1 and (68 1, respectively (with h = 1). The Luttinger parameter for the liquid is taken to be «: = n/mc = 20, while 
M = m. 
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VI. FINITE TEMPERATURE DYNAMICS OF DEPLETON. BACKSCATTERING OF THERMAL 

PHONONS AND VISCOUS FRICTION 



At a finite temperature, even in the absence of the external force, the collective coordinates fluctuate around their 
equilibrium position Aj?^ , found from the condition d^H = 0. Assuming these fluctuations are small, one may linearize 



Eq. (46 1 near its equilibrium point. Being transformed to the frequency representation, such linearized matrix equation 



of motion takes the form 



Y 



A(c.) - en 



(69) 



where V is the Hessian matrix of the second derivatives of the impurity Hamiltonian at fixed momentum P. Its matrix 
elements are 
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Ha-A+ Ha_a- 



(70) 



A=A(0) 



Solution of equation ( 69 ) takes the form 



Aiu;) = [i-'-^r\ re(c.) = re(a.) + yr2e(^) 



(71) 



where we consider it as a perturbative sequence in frequency, to avoid spurious instabilities mentioned in Section [V] 
Substituting this solution into the right hand side of Eq. ( 45 1 and averaging it over the Gaussian noise ( 43 ) , one finds 
for the momentum loss rate 
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cLu 
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w^Tr 
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(72) 



where we kept only leading order in frequency. The matrix-valued function F(cj) is odd in frequency, selecting only odd 
powers of uj from the expression it is multiplied by. Equation ( 72 ) provides the expression for the viscous friction force 
acting on the impurity from the normal component of the liquid. It can be identified with the Raman two-phonon 
scattering mechanism [TUl HH ESI HI] . 



Substituting the explicit form of the Ohmic noise correlator (41) and performing the frequency integration, one 
finds for the friction force 
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(73) 



The T'^ dependence of the friction force was reported in [T3] and is a direct consequence of the phase space for 
two-phonon scattering. It can be also viewed as a one-dimensional version of the Khalatnikov-Landau result [TU] 
for the viscosity of liquid helium. It should be noted that the friction force ( 73 ) is proportional to the off-diagonal 

matrix element F-i which represents backscattering of phonons by the depleton. It depends on the momentum 

P and, consequently velocity V of the mobile impurity in addition to the parameters of the impurity-background 
interactions. Below we derive the general expression for the backscattering amplitude and relate it to the equilibrium 
dispersion E{P, n) of the depleton. 



A. Backscattering amplitude 



The detailed information about microscopic impurity-liquid interactions is contained in the off-diagonal matrix 

element which represents the effective vertex for two-phonon scattering. In the spirit of our phenomenological 

approach it may be expressed, as explained in Appendix [Ej through partial derivatives of the equilibrium values 
MV,n) andiVo(l^,n), 
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(74) 



Treating collective variables $0 = ^Q{P,n) and A^o = No{P,n) as function of momentum rather than velocity and 
calculating the corresponding derivatives as explained in Appendix |E] we find an alternative representation 



(75) 
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FIG. 7. V — and $ = (i.e., P — 0) dissipation rate as a function of the impurity coupling G for n ■ 
line). The dashed line corresponds to the strong coupling G ^ c limit, Eq. (761. 



10 and M = m (solid 



The functions $o A^o ^tre obtained from the equilibrium dispersion E(P,n) as explained in Section IIIB (see, e.g. 

Eq. (21 1). This relations can be used to obtain an equivalent representation for the backscattering amplitude 

given by Eq. (Ell I in Appendix pj 



Two remarks are in order. First, it should be noticed that the backscattering amplitude given by Eqs. (74), (75) 
depends on the velocity V or momentum P of the depleton. While the momentum dependence is unambiguous, 
there can be different regimes corresponding to the same velocity characterized by different values of the superfluid 
phase Second, it can be shown by calculating dispersion E{P,n) by Bethe Ansatz method that the off-diagonal 

matrix element T-i vanishes identically for integrable mobile impurities, such as lowest excitation branch in Lieb- 

Liniger [JS] or Yang-Gaudin [TS] models and the soliton of the Gross-Pitaevskii equation (see Appendix [b| . We 
postpone the details of these exact calculations to another publication. Below we use the semiclassical approach valid 
for the weakly interacting bosonic background which provides results for backscattering amplitude in a wide range of 
parameters relevant to experiments. 



1. Strong coupling regime G/c 3> 1 



In the strong coupling regime the leading terms in the the number of depleted particles, the second equation in (C9) 
is a constant, A^o = 2k ^ 1 and the momentum is dominated by the supercurrent and P = n^. These observation 
simplify greatly the expression ( 75 ) for backscattering amplitude and one gets 



1 



1 



M 
Km 



(76) 



where we have used the fact that dn/dn ~ n/2n. The expression (76) is momentum or velocity independent constant, 
which equals for a not too heavy impurity M ^ nm. 



2. Weak coupling regime G/c <^ 1 for slow impurity 



In the weak coupling regime the backscattering amplitude varies strongly with the velocity of the depleton. Here 
we present results for slow impurity which can be found by evaluating r_| for V — 0. 

There are two distinct physical regimes corresponding to a slow impurity: one, corresponding to the solution with 
$0 ~ and momentum P ~ MV ~ dominated by the bare impurity and one corresponding to $o — and 
momentum P ~ irn dominated by the background supercurrent. For <I> « we have 




$o(V,n) = 2- 
c 




(77) 
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while near $ = tt. we find 



V 



1 



(78) 



Tliese expressions are correct up to terms 0{V'^) for the number of particles and 0{V^) for the phase. 

In the case $o — we use Eqs. (77 1 to calculate For weak coupling to the leading order in both g/c = 1/k 

and G/c we have, 



1 



G/c< 1. 



(79) 



This result agrees with Eq.(5) in Ref. |14] and indicates that friction vanishes near integrable point M = m and 
G = g, at least for small velocities near $ = 0. The results for arbitrary coupling G/c are presented in Figs, [t] 
In the case $o — ^r, we use Eqs. (78) and calculate effective mass M* = M — niN+ndv^ = M — 2Km — 2nm{\ + G / c). 



Using Eq. ( 74 ) we find for arbitrary coupling 



1 {M/Km){l + 2G/c)-2G/c 
^ M/Km-2G/c- 4 ' 



(80) 



It comes as no surprise that Eq.(80) vanishes in the limit M, G — 0, where one recovers the dark soliton result 



from Appendix |B] Indeed, there it is mentioned that such soliton excitation is transparent for phonons due to its 
integrability and, consequently, F_| = for all velocities. In the limit of weak coupling Eq.(80) becomes 
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(81) 



At the Yang-Gaudin integrable point (|T5J, see also 



25j ) of the quantum system defined as M = m and G = g, 
the corresponding limit of the second equation in (81) is proportional g/c = 1/k which is a small parameter in the 



semiclassical limit and therefore our calculation goes beyond the accuracy of the approximations adopted above. 



B. Bloch oscillations in the presence of viscous friction 



Imagine the impurity is dragged by the weak external force F. There is a velocity V < Vc such that F + Ffy{V) = 
and the system reaches a steady state with constant drift velocity Vb and no oscillations, Vb = 0. If the applied force 



is very small, the drift velocity can be found from the small velocity limit of Eq. ( 73 ) in the linear response form, 
axF^ where the finite temperature mobility is given by 
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|F+_(F = 0)| 



(82) 



in accordance with [131 El] ■ Increasing the external force, leads to a slightly nonlinear dependence of the drift velocity 
on F (due to the non-linear velocity dependence of ^fr(y), cf. Eq. (73)), until the the maximal possible velocity Vc is 
reached for F = i^min = —Fi^{Vc). Beyond this critical force no steady state solution can be found and the impurity 
performs Bloch oscillations along with the drift. In such a nonlinear regime the amplitude, period of the oscillations, 
drift velocity, as well as the momentum-dependent viscous backscattering amplitude, Eq. (75), are controlled by the 
equilibrium dispersion relation. 

We illustrate the above scenario using the model of a strongly coupled impurity. In this case the critical velocity is 
small and backscattering amplitude is velocity-independent. The critical force is then found from the linear response 
as Fniin = Vc/(tt- We can use the relation V — VcSin^ + K$/2n (i.e., the second of Eqs. (61)) and the fact that 
P = n$ to write down the equation of motion for the superfluid phase. 
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where A — F/F^i^. Introducing dimensionless time variable s — {F^ 
of generality, $(0) = 0, the solution of Eq. (83 1 is found to be 
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(84) 
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FIG. 8. Drift velocity of a strongly coupled impurity in a weakly interacting bose liquid as a function of temperature T 
for various forces as predicted by Eq. (86 1. F is given in units of -FmaxK/c, while the Luttinger parameter is taken to be 



K 



■ 52 so that, to a good approximation, -F'min(T - 



max ' 



For A < 1 it describes velocity approaching its limiting value Voo = Ksin$(oo) — XVc- For A > 1 Eq. (84) describes 
a periodic function with period 

t 277 27rn(l + n/2aTn'^) 

The corresponding drift velocity for F > i^min can be found by averaging the momentum relation n<i> = F — V/ctt 
over a single period with the result 



Fo^vi-^-A(?2SIZiy (80) 



For Fmin <^ F ^max, One may expand Eq. (86) to obtain a small temperature correction 



which gives the previous result Eq. (59) in the limit of small T, or more specifically, small Fmin- Interestingly, Eq. (87) 
predicts an increase in the drift with increasing T if -F/-Finax ^ Vc/c. This suggests that for a small enough force, 
a significant drop in the drift velocity may occur as the system is cooled below the critical temperature, when the 
impurity enters the regime of Bloch oscillations, see Fig. |8] 

The above analysis demonstrates, inter alia, a non-monotonic dependence of the drift velocity on the parameter 
F/Fjnin as it is increases past 1 and enters the regime of Bloch oscillations. This may occur either by fixing the 
temperature and increasing the force, or by fixing the external force and cooling the system. In contrast to the 



vanishing amplitude Vb near i^max, indicated by the second equation in (64), it is rather the divergent period, Eq. (85) 



for F approaching Fj^in from above which leads to the disappearance of Bloch oscillations 



VII. DISCUSSION OF THE RESULTS 



The coupling with phonons governed by the universal term, Eq. (29), results in transfer of energy and momentum 
between the depleton and the background. At zero temperature it takes the form of radiation of phonons by accelerated 
impurity similar to the radiative damping in classical electrodynamics |40j and leads to the finite mobility of the 
depleton. For sufficiently weak forces the mobility, defined as the ratio of the drift velocity to the applied force, see 
Eq. (56 1 can be expressed via equilibrium dispersion, Eq. (57) reminiscent of the Kubo linear response theory. The 
drift, however, is superimposed with the essentially non-linear Bloch oscillations of the velocity. The amplitude of 
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FIG. 9. Schematic force vs. temperature diagram of impurity motion in a ID quantum liquid. In the region F < Fmin (hght 
gray) the impurity drifts with mobility or oc T'^ for r < mc . For T ^ mc , the mobility scales or oc T"^ (dashed are the 
small and large temperature asymptotics). For -Fmin < F < -Fmax (dark gray) one has Bloch oscillations + drift. Above Fmax, 
our theory is inapplicable and we expect some kind of incoherent acceleration, possibly corresponding to a supersonic impurity 
which has escaped its self-induced depletion cloud. 



the oscillations is shown to vanish as the external force attains the upper critical force -Fmax reflecting the limit of 
validity for the description based on the scale separation between the healing length and the phonon wavelength. 

At finite temperature the thermal phonons present in the system are scattered by the depleton leading to the viscous 
friction force, Eq. (73). This in turn leads to the appearance of the lower critical force Fmin(r) ~ which sets the 



lower limit for the external force driving the Bloch oscillation. In contrast to the situation in the vicinity of Fmax, the 
approach to Fmin makes the Bloch oscillation disappear through the divergence of the period, Eq. ( [SS] ). Below Faun 
the system enters the regime of non-oscillatory drift characterized by the velocity Vu < Vc-, where all the momentum 
provided by the external force is dissipated into phononic bath. 

The viscous force contains valuable information about the fine details of the interactions between particles. In this 
work we have confirmed and extended the earlier observation [14j, that the backscattering amplitude, given in terms 



of equilibrium properties by Eqs. ( 74 1 , ( 75 ) vanishes if the depleton is an elementary excitation of an integrable model. 
This includes the dark soliton excitation of Lieb-Liniger model [35| as well as spinous in bosonic Yang-Gaudin model 
[15j . The microscopic mechanism responsible for the absence of dissipation is due to destructive quantum interference 
between various two-phonon processes. It can be traced back to the absence of three-body processes lying at heart 
of integrability in one dimension. In contrast, the radiative processes due to the presence of external force are always 
present in the dynamics of the depleton. This is because the external potential is not, in general, compatible with the 
integrability. 

Various dynamical regimes of depleton are summarized in diagram Fig. [9j At sufficiently low temperatures there 
is a wide parametric window i^min < F < i^max for the external force in which the Bloch oscillations can be observed 
experimentally. At temperatures higher than chemical potential the mobility of the particle becomes inversely pro- 
portional to temperature [43l |44] and moderates the growth of F^i^ (see Fig. |9]). The above range of forces can be 
increased further by exploring the dependence of Fmin and -Fmax on the interaction parameters. In particular, at 
integrable point -Fmin vanishes for any temperature. 

Our approach to the dynamics of depleton is essentially classical, the quantum mechanics enters only via parameters 



of the effective action. The equations of motion, Eqs. (46) neglect therefore the quantum and thermal fluctuations 



of the collective variables of the depleton. These fluctuations can be taken into account by either simulating the 
Langevin equation in the presence of equilibrium noises or by writing appropriate Fokker-Planck equation for the 
distribution function. One of the most important consequence of the fluctuations is expected to be the smearing the 
boundaries between dynamical regimes in Fig. [9| We leave this as well as the question about the role of quantum 
fluctuation for further investigation. Another important extension of the present work would be the studies of the 
finite size effects due to the trap geometry relevant for experiments with ultracold atoms. 
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Appendix A: Derivation of dissipative action 



Gaussian integration of the interaction term Eq. ( 38 ) using the following phononic propagators 



- »(xci(^, t)xUx', t')) = D^(a; - x' , t-t'), -i{x^ix, t)xUx' , t')) = D^(.x -x',t~ t') 
-il^xA^. t)x\{x',t')) = D«(a; ~x',t- t'), -i[x^{x, t)x\{x' , t')) = 
leads to the quadratic nonlocal action 

5eff dm' Alit) dt [A«(t - t') - A^{t - t')] [Aq(t') + V-iAci(i')^q(i') 



(Al) 



dtdt' 



AUt) + x(t)Al{t)y-'] d^A^it t') k(t') + y-^A,,{t')x^{t') 



(A2) 



where A^'^'^(i) — D'^'^'^ {Vt^t) are phonon propagators restricted to the impurity trajectory. 

Inverting the matrix in Eq. ( 33 ) and taking appropriate analytic structure in the complex w plane we obtain Fourier 



components of the retarded and advanced propagators, 

t 



D"{q,Lo)= D^{q,uj) = dxdte-"^''+"^'D^{x,t) 



This leads to 
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The noise terms, i.e. the second line of Eq. ( A2) are controlled by the Keldysh component and its derivatives. 

m we 

(A6) 



Assuming thermal equilibrium of phonons in the laboratory frame and using Fluctuation-Dissipation Theorem we 
have 



D^iq, u) = coth (^) (D«(g, u) - D^(q, u)) 



Using Eq. (A4) we find the Fourier component of the Keldysh propagator restricted to the classical trajectory 
dq 



2^ ( ^^^^ ) ( D""!?, ^ + qV)~ D^(g, u + qV) 
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Substituting its Fourier transform it into Eq. ( A2 ) and taking double time derivative leads to the second term in 
Eq. (|40l). 



Appendix B: Grey Solitons 

Dynamic properties of one-dimensional bosons of mass m weakly interacting via repulsive short range potential 
proportional to coupling constant g can be studied within the Gross-Pitaevskii description using the following La- 
grangian 



L = dx [ iipdttp 



f I^I' + mIV'I' 



(Bl) 
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Here is the quasi-condensate wavefunction corresponding to the asymptotic density \ip(±cx),t)\'^ = n and vanishing 
supercurrent at infinity. The condition of weak interactions is mg/n <C 1. Minimizing the action defined by the 
Lagrangian Eq. iBll leads to the Gross-Pitaevskii equation, 



2m 



(B2) 



3^ related to the sound 



Substituting a constant solution ipQ — ^/n one obtains the chemical potential fi — gn ~ 
velocity c = yj^ffnfm. In addition to uniform solution, Gross-Pitaevskii equation (CI I admits a one-parameter family 
of solutions 



i\)s(,x — Vt; $s) = vn cos 



tanh 



x-Vt 
I 



(B3) 



known as grey solitons [28l |3T] . They can be visualized as a dip moving with velocity V and having a core size 
I = l/m^y (c^ — V^). The solution, Eq. (B3) is characterized by the total phase drop and number of expelled particles 
related to the velocity as 



= 2arccos^, Ns{V,n) ^ j dx {n ^ \Mx)f) ^-^c^-V^. 



Momentum and energy of the soliton are given by [28l HH |49] 

Ps = - mN^V , 
which allows to define the Lagrangian 



i?. = ^cnsin3(ci>,/2)="^iV3 
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Using soliton Lagrangian (B6| and putting M = in Eqs. (17 1 and (18 1 one sees immediately that the variables in 
(B4| coincide with the collective variables $, TV. Therefore the grey soliton can be viewed as a model for a massless 
impurity consisting of depletion cloud only. 



We invert Eqs. (B4) which yields 
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Using Eq. (20) together with Eqs. (B4| and (B6| yields the internal energy, Eq. (23) of the soliton, 
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The matrix of second derivatives at the equilibrium solution reads 
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The backscattering amplitude r_| , calculated with the help of Eq. ( |E1| , vanishes identically due to the integrability 

[50] of the Gross-Pitaevskii equation (B2|. 

As it was shown in (44i A6\ a weak cubic nonlinearity —{a/6)\ip\^ in the Lagrangian (Bl) breaks the integrability 
of the model. Cubic terms describe three-body interactions which arise from virtual transitions to higher transverse 
states of tightly confined one-dimensional liquid 02 ■ Here we extent the results in Ref . 02 and calculate the amplitude 
of the corresponding dissipation processes for the whole range of soliton velocities. 

The corresponding correction to the Lagrangian, Eq. ( B6 ) of the soliton can be calculated to the leading order by 
evaluating it with the unperturbed solution. 
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FIG. 10. Superfluid density in the presence of an impurity moving with velocity V/c — 0.6 and coupling G/c = 0.1 (solid hne). 
The dashed lines corresponds to the double sohton solution Eq.(B3l used to construct Eq.(C2|. 



Here we have used the expression Eq. ( B4 ) for the number of expeUed particles N to calculate the correction to the 



energy relying on the theorem of small increments. The corresponding change in the matrix ( B9 ) of second derivatives 

^H--^^^x(0;) (BID 
can be taken perturbatively in the calculation of the off-diagonal matrix element JF-i We have 



r T-iJH (T-1)' r 



Substituting Eqs. ( |B11| , (jSlj), yields 
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Here we have used the fact that the matrix f is diagonal in the leading order in a and det£T^T) = 1. The determinant 



of the matrix (B9) is det H — (rnc^ sin (f>/2) . Substituting it into Eq. (B13) and using (B4) leads to 
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in agreement with the results of Ref . [43, 



Appendix C: Impurity in a weakly interacting liquid 



To model the impurity coupled to a weakly interacting superfluid at T = 0, the Gross-Pitaevskii equation, Eq. (B2) 
is modified in the presence of a delta function potential moving with constant velocity V 



' 2m ' 



(CI) 



For G > 0, the soliton solution, Eq. nBSl still satisfies Eq.(Cl ) except at the location of the impurity. Thus, one may 



construct a solution to Eq.(Cl) by matching two soliton solutions, Eq. (B3) at the location of the impurity, as shown 



in Fig. 10 The proper solutions of Eq, (CI ) can thus be written 
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where y = x — Vt and velocity is always related as V/c = cos($s/2) to the phase <J>s parameterizing the soliton 
configuration, Eq. (B3). The solutions 'ijj± to the right (left) of the impurity satisfy the boundary conditions: ip+(0) = 
'0-(O), ^^+(0) — V'-(O) = 2toGV'(0). Using (C2|, the boundary conditions give rise to the following two equations for 
$0 and z = tanh(a;o/0 
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Equations (C3|, (C4) permit a solution only for V < Vc{G/c) where Vc is some critical velocity that depends only 



upon the parameter G/c [T71 ST] . This can be seen by considering the right and left hand sides of the second equation 
( C4 ) . While the left hand side is bounded by the maximum at Zmax = 1/ the right hand side grows quadratically 
with z and therefore solution exist only for a limited range of $s, which leads to the above-mentioned limitation on 
velocity. 



For this reason we choose to parameterize the solution, Eq.(C2|, by the total phase drop across the impurity, 
$ = $s — which happens to permit a solution for any $. Thus, upon solving Eqs. (C3), ( |C4| one finds the 
relations z = z{^,G/c) and $s = ^s{^,G/c). It can easily be seen from Eqs. (C3), (C4| that these functions are 
periodic in The number of expelled particles and momentum can be calculated and expressed through the phase 
$ as 

dx{n - IV'P) = 2k sin %(1 - z); P = n<I> + (M - mN)V, 
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The energy may also be calculated and expressed in terms of <f> as 
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Alternatively, we may solve for (f> = by inverting the second of Eqs. ( |C5[ ). Substituting it into the energy 

function, Eq. (C6) one obtains the dispersion relation E{P) plotted in Fig. [I] for the impurity in a weakly interacting 
bose liquid. 

In the weak coupling regime G/c 1 the critical velocity can be obtained by neglecting the z-dependence of the 
r.h.s. Eq. ( |C4[ ). This is justified a posteriori as at the solution 4>s ^ 1- Expanding the trigonometric functions and 
using Zmax = I/a/S in the l.h.s one obtains ($8/2)"^ = (3\/3/2) G/c, which justifies our approximation. We thus have 
the critical velocity 

-^-li — T- (C7) 



V, 



Solving Eqs. (C3), (C4) for arbitrary velocity or momentum is cumbersome and we resort to numerical methods. 



In the strong coupling limit G ^ c we determine the dependence of z and <I>s on <& to order c/G. This is done by 



writing <I>s — tt 



;$s^^ and z 



and finding coefficients ^s^'' and z^^^ from Eqs. (C3), (C4|. We have 



z($) ~ ^cos'^; 



Using Eqs. (C8|, (C5) one has 
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The first equation has the Josephson form with the critical velocity Vc = c^/2G in the strong coupling limit. It is 
also clear from the second equation that in the leading approximation is a constant A^ ~ 2k. This results in the 
Josephson form for the energy, E q. pi] ). 

At V = the equations (C3), (|C4| simplify considerably by putting $s — t^- There are two solutions. 
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2c 



0, 



(CIO) 



The zq root corresponds to $ = and describes a background only slightly perturbed by the stationary impurity. The 
z^ root corresponds to the dark soliton solution with <I> = tt, which persists for G > because the density vanishes at 
the impurity location i.e., there is no additional energy cost to put the impurity in the center of a dark soliton. 
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Appendix D: Solution of equation of motion for a strongly coupled impurity 



We write the sinusoidal forms for the velocity and the phase of depleton 

V{t) ^VB + VBam{ujBt + 6v); ^{t) = a;B(l + A cos(wBt + 5$)), 



(Dl) 



depending on a priori unknown parameters Vdj ^Bj (^Bi (^Vi <^*j ^ and substitute them into Eq. (61 1. One first arrives 
at the following relations between the time independent components. 
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Solving these equations one obtains Vb and wb hi Eq. ( 64 1 . In the limit F ^ i^max — 2nmc^ we recover the linear 

der de 



dependence Vb/c = F/F^i^^. The leading order deviation from linearity is given by 

F / F'^ 



F„- 



1 



(D3) 



Substituting the ansatz (Dl ) into the first of Eqs. (61 1 gives a relation between the time dependent components, 
KwlVb 



A 



noJB 



cos(wBt + 6^) + Vb{M - mN)ujB 



cos(^Bt + 5v) + 2 ^ {c^~V^Y + '^^^J ^ 

(D4) 

where we neglected terms 0{A^), 0{V^) and 0{AVb) to keep the calculation to first order in Vc/c, as both Vb and 
A will be seen to scale with Vc/c. The second term in brackets in Eq. (|D4[) is simplified employing the formula 



cos X + a sin 



X = —\/l + oP cos{x + 7r/2 + arctan(l/a)). 



(D5) 



In order to cancel the term cx cos(ajB^ + <^*), we require a definite relation between the phases and amplitudes. The 
constraint for the phase is Sy + 7r/2 + arctan(l/a) = (5$. From Eq. (D4) we have 



1 kujb 



2 M - mN [c"^ - V^Y M -mN {l-V^/c?Y c 
The equality of amplitude implies a relation between A and Vb, namely 



A^ 



Vb M - mN 



1 + Vg/c2 



Finally, we substitute the ansatz into the second of Eqs. (61 ) to obtain 

"^B 1 / , , r N , 



AVn 



cos{ujBt + S^) + 



VT 



aiii{ujBt + (5$) — -Vc sin(cjBi + $0), 



(D6) 



(D7) 



(D8) 



where we used sin(a;Bi + iy) = —-j== (sin(wB^ + <^$) + acos(a;Bi + (^$)) and $0 is some initial phase of $ which 



comes from integrating the second of Eqs. (Dl). Using Eq. (D5) again finally gives 

\-V^lc^ 



VB^v^^m^^v 



v/TTv^T^ ^(1 + v%lc^f + {M^y (1 - V^/c^f V^/c 



(D9) 



where Vb is given by Eq. (64). As expected, Vb is an even function of F since Vb is odd. For small F F„ 
Eq. (IDQI) gives Eq. $^. 



Appendix E: Calculation of the backscattering amplitude r4 



Using the transformation Eq. (35) we have f = T^H which leads to 

r+- == ■ 1 1 { kHmn 

det H V K 



(El) 
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The matrix of second derivatives 



H 



(E2) 



dependence of the velocity y on $ and N at constant momentum P, one can rewrite Eq. (E2) as 

,2 



is calculated at equilibrium values of <i> and N. Differentiating Eq. (11) and taking into account Eq. (10 1 for the 

(E3) 



H = 



1 



M - mN 



n" —nmV 
-nrnV ■w?V'^ 



1 



detQ 



The last term in this equation represents the matrix of second derivatives of H^{'^,N). Owing to properties of 
Legendre transformation we have expressed it as the inverse of the Hessian matrix, 



(E4) 



of the thermodynamical potential Hereafter we drop the primes over $7, j and fj, for clarity. In writing 

Eq. ( E4 1 we used Eq. ([5| to express double derivatives as derivatives of equilibrium values <1> and N with respect to 
underlying values of the supercurrent and chemical potential. 

dfi/dV, n 



Using relations mV 
detH = 



1 

detQ 



1 

M - mN 



dj/dV in the first term of Eq. (E3) simplifies considerably the determinant 

M* 



[m^V^Nf, + - mnV{Nj + $^)] 



1 



detQM - mN ' 



(E5) 



by using the effective mass Eq. (19 1 of the equilibrated impurity. We use this fact and invert the matrix in Eq. (E3) 



by a standard procedure and find for the off-diagonal matrix element 

1 



kM* 



(M - mN) {k^^j - N^) + (1 - V'^/c^) ($jiV^ 



(E6) 



In the heavy particle limit M ~ M* ^ m only the first term in the square bracket survives. In this limit can 



be identified with backscattering amplitude of a static impurity F-i 

To obtain this result we inverted the relations in Eq. (16), 



3 - 



1 



m (c^ 



'—mc^/n mV 
y2) V -V n 



dv 

dn 



^N^, = - + mdnN) /mc. 

(E7) 



and used the identity n9„$ — mVdnN = Vdv^ — {mc^ /n)dvN . The latter is obtained from equality of the mixed 
derivatives of the Lagrangian L{V,n), obtained by differentiating Eqs. (17 1, (18) by n and V respectively. The 
second term in the square brackets in Eq. ( E6 ) is transformed by applying ( E7 ) . Combining it with the term, 
(Af - mN)/M* yields Eq. 



{d/dn)^ = (d/dn) 



To deal with equilibrium functions $(P, n), N{P,n) obtained for a given momentum P rather than velocity V we 
use the fact that dv — M*dp and the following relation for the derivatives with respect to the density 

M* {dV/dn)p dp . (E8) 

Differentiating Eqs. (17), (18) with respect to P and n at constant momentum leads to 

1 / m(? jn mV 



dpN 
d„N 



m (c^ 



mc^/n mV 



m (c2 -V^)\ V 



(M* - M + mN)/M* 
dnV + <i>/M* 

- (Af - mN)dnV 
$9„F + {aa - aN) 



(E9) 
(ElO) 



where ad = d^E. Substituting Eq. (E8) into Eq (74 1, using the effective mass, Eq. (19) together with Eq. (E9) we 
arrive at Eq. ( 75 ) . Alternatively one may wish to express directly in terms of derivatives of the dispersion law 



E{P, n) with the result, 



1 - 



M - mN 



H ^(a<j - aN) 



{M - mNY - m^K^^^ 
mM* 



(Ell) 



In the limit of small velocities terms proportional to $ can be neglected and this expression becomes is identical to 
Eq. (3) in Ref.jM] with N =\-bi. 
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